NEW  YORK  UNIVERSITY 

IMM-NYU  169 

INSTITUTE  OF  MATHEMATICAL  SCIENCES 

LIBRARY 
13  V/avarly  Place,  New  York  3,  NY. 


THE  DOCK  PROBLEM 


Bt 


K.  O.  FRIEDRICHS  and  HANS  LEWY 


o— 
'-9 


'     —  Prepared  under 

^  (J  Navy  Contract  N6ori-201  Task  Order  No.  1 

M  By  the 

Institute  for  Mathematics  and  Mechanics 


New  York  University 

October,  1947 


The  Dock  Problem 
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Suppose  one  half -plane  of  an  infinite  water  siirface 
is  covered  with  a  rigid  plate,  the  "dock,"  How  does 
this  dock  influence  waves  standing  or  traveling  perpen- 
dicularly to  the  edge  of  the  dock  on  the  free  water  sur- 
face? This  problem  may  be  considered  as  a  special  case 
of  the  problem  of  waves  on  a  sloping  beach,  viz,  for  a 
beach  with  a  slope  angle  of  180  ,  see  the  introductory 
report  [4],  For  beach  angles  which  are  an  integer  fraction 
of  a  right  angle,  ^^  ~  "^  *    ^^^   problem  was  solved  by 
Miche  [1],  Stoker,  and  Lewy  [2,3];   for  beach  angles  of 
the  form  co  =  £2  with  relatively  prime  integers  p  and  q, 
the  problem  was  solved  by  Lewy  in  case  p  is  odd.  The 
present  problem  corresponds  to  the  case  q  =  1,  p  =  2. 
Although  it  may  be  possible  to  make  ^  tend  to  the  limit  1 
and  follow  throiigh  this  passage  to  the  limit  in  the  explicit 
expressions  for  the  waves,  we  attack  the  problem  differ- 
ently, using  a  method  related  to  the  Laplace  transforma- 
tion.  It  turns  out  that  the  problem  permits  a  simple 
explicit  solution  which  offers  the  possibility  of  discussing 
the  nature  of  the  solution  and  of  determining  it  numerically, 
We  shall  in  particular  discuss  the  behavior  of  the  wave 
motion  near  the  edge  of  the  dock,  i.e.  near  the  line  along 
which  the  water  surface  and  the  dock  meet.   In  addition, 
we  shall  show  graphically  the  amplitudes  of  the  waves  on 
the  free  surface  and  of  the  pressvire  under  the  dock. 

Aside  from  its  interest  as  a  flow  problem,  the  dock 
problem  is  of  significance  as  a  new  example  of  a  potential 
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problera  with  different  linear  boundary  conditions  on 
different  parts  of  the  boxmdary,"" 

'    Since  we  are  interested  in  two-dimensional  flow, 
all  quantities  are  assumed  constant  along  any  parallel 
to  the  edge.  We  place  an  (x,y) -coordinate  system  such 
that  the  water  fills  the  lower  half -plane  y  <  0,  the  free 
water  surface  is  the  positive  x-axis,  y  =  0,  x  >  0,  and 
the  dock  is  the  negative  x-axis,  y  =  0,  x  <  0,  We  set 
X  +  iy  =  z  =  r  e   •  Further  we  choose  the  xmits  of  length 
and  time  such  that  the  frequency  of  the  oscillations  of 
the  water  surface  and  the  acceleration  of  gravity  have 
the  value  one.  The  water  flow  is  then  described  by  a 
complex  potential  function 

\{z)   cos  t 

from  which  the  velocity  components  u  and  v  are  derived: 

u  -  iv  =  ^^^^  cos  t  • 

The  boundary  condition  that  the  pressure  be  constant 
on  the  free  surface  leads,  by  the  well  known  process  of 
linearization,  see  e.g.  Stoker  [3, p. 6],  to  the  condition 

^      ^^4  "X  +  iX)  =  0,   y  =  0,  X  >  0  , 

while  the  condition  that  the  water  slides  under  the  dock 
siirface  leads  to  the  condition 


"A  solution  of  the  problem  of  two-  and  three-dimensional 
waves  over  a  channel  of  finite  depth  with  a  dock  has 
been  annoxjnced  by  A.  E.  Helns  [5]. 
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II  Im(^  X)    =  0   ,        7  =  0,  x<0. 

At  infinity,  i.e.  as  [z]  — >  oo,  we  further  require 

III  X(z)  -^  x^ 

i.e.      "^(z)e^^  ->  "X.^   ,   0.,  <  Q  <  0  , 

"X.(z)  ->  0  ,     -1T  <  e  <  e.,„  , 


for  any  angles  9...., 9,.,_„,  with  -tt  <  0.,  <  ©„_„_  <  0  and  with  an 
appropriate  complex  constant   A   ^0. 

These  two  conditions,  ?/hich  are  evidently  not  contra- 
dictory, assure  that,  at  large  distances  from  the  dock, 
the  water  executes  sinusoidal  motions  near  the  surface 
and  quiets  down  in  other  directions. 

At  the  edge  we  permit  a  logarithmic  singularity,  i.e. 
we  require  that,  with  an  appropriate  constant  k, 

IV      XCz)  -  k  log  z   is  bounded  at   z  =  0  . 

For  the  significance  of  this  logarithmic  singularity  see 
Stoker  [3,  p.  11] . 

We  shall  show  two  solutions  of  the  problem:   one, 
■p 
7^  (z),  bounded  at  the  edge  z  =  0,  the  other,  given  by 

(1)  A^(2)  =-A  X'^C-) 

and  possessing  a  logarithmic  singularity  at  the  edge. 
Specifically,  the  beliavior  of  the  solutions  at  z  =  oo  is 
given  by 
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iTT  ^ 

(2)  72  X?(z)-  e"^  e"^^  ->  0,n/2X?(z)-  le^  e"^^  ->  0,  z  ->  oo 

if  z  approaches  infinity  within  a  sector  which  excludes 
the  negative  real  axis,  -rr  <  0_„  <  ©  <  0,  while  if  z 
approaches  infinity  in  a  sector  which  excludes  the  positive 
real  axis,  -tt  <  ©  <  ©„„  <  0, 

(3)  X^(z)  -*  0,   X^(z)  ->  0,   z->oo. 

Both  statements,  (2)  and  (3),  are  in  agreement  in  the 
overlapping  parts  of  the  sectors.   The  behavior  at  1±ie 
origin  is  given  by 

(4)  ^iz)    ->  I   as   z  ->  0  , 


(5)    X^(z)  -  -k   ^°S  (-Z)  ->  -  ^  as  z  -i^  0  , 


C  being  the  Euler  constant. 
The  relation 

TT 

(6)  X^(2)e^^P/x^(0)  =  v^  e  ^ 

is  particularly  remarkable,  since  it  is  in  agreement 
with  the  value  "Y^  \/p/q  found  for  the  same  ratio  for  beaches 
of  slope  angle  PTr/2q,  see  [2, p.  744],  Y^  being  an  appropriate 
number  of  absolute  value  one.  Note  that  in  o\ar  case, 

p=2,  q=l,  pis  not  odd. 

R         S 
The  solutions   X  (z),   X^^^  a^®  given  by  the 

integrals 
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(8) 


X^(^)  =-^/  •^^  C1^(C)  dC  . 


in  which  (P    is  any  path  running  in  the  lower  half-plane 

■n  <  0  from  4^  =  0  to  C  =  oo  entering  C  =  c»  in  a 
direction  along  which   Jz  has  a  negative  real  part  and 
skirting  the  point   4^  =  -1  so  as  to  leave  it  to  the 
right.  Such  a  path  (r  will  then  be  called  "admissible" 
for  the  value  e* 

The  fimction  V  (^ )  is  defined  in  the  lower  half- 
plane,  ^  <  0>  (^  /^  0),  except  for  the  point  C  =   -1# 
at  which  it  has  a  pole.  It  is  given  by  the  two  equivalent 
representations 

_  1  ^1 

(9)o  V  (C)  =  ^ "  ^(1  -  ir)"^(i  + 14:)"  ^ 


if  1      ^^ 

exp  -  J     log    (iC^)-     -5-—.  , 


^  1 


n-xoo  ^^^ 

exp  -  J  log   (1^^)  -2—-  . 


^-lOD 


dC. 


Note  that  the  fimction  /    log  (i^„)  —n — ^  is  regular 
in  the  lower  half -plane.  The  equivalence  of  the  two 


-6- 


representations  follows  from  the  relation 


(10) 


4*  +  l 


.00 


■K- 


Prom  the  relation 

-1  .y  «-l00 


f    log  (i^  )^Jt-  =  r       log  (ri)4^ 


+  1 

obtained  by  substituting  -  i  for  ^^  and  from  (10),  we 
conclude 

log  (ir„  ^ 


^^^)         /   ^°S  (l?,)-2^  =  ^ 


and  we  see  that  (^  +  i)  ^(^)  is  regular  at   ^  =  -i 
and  assvimes  there  the  value 

,-i      ^^ 
(12)  (4+  i)  ^{C)\        =  1  e"^  . 

The  path  of  integration  ^  is  any  curve  coming  from 
the  origin,  skirting  the  point   ^  =  -i  leaving  it  to  the 
right,  and  then  going  to  infinity  in  a  direction  in  which 
the  real  part  of   ^z  is  negative.  Thus  from  the  ex- 
pressions (9)^  and  (9)QQit  is  clear  the   ^(4^)  behaves  at 
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5=0  and  at   ^  =  oo  in  such  a  way  that  the  two  inte- 
grals (7)  and  (8)  represent  analytic  ftmctions  of  z  which 
are  regular  in  the  lower  half -plane,  y  <  0,  (z  ^  0), 

In  order  to  verify  that  the  boiindary  conditions  I 
and  II  are  satisfied  we  use  different  representations  of 
the  function   V(C)«  Prom  the  relations 


,c 


(13)^ 


i  /  -^  =  ^   log  ^ij|  ,    except  for  4  =  0,  1^1  >  1, 
o  C-x- 


vOO 


4 


we  find 


_  1       _  3        „l  ^y 

(14)^   V(C)  =  ^"  ^(l+<^)"  ^  exp  i  /   log  ^-^  -5-ii:  , 

(15)^   V'(^)  =  r  '^(l-iC)"^(l+C^)"  ^  exp  i  /   log(-C,)"^  -2-^  * 


valid  for  t1  <  0  except  for  <f  =  0,  y|  <  -1, 

(14)^^  V^(C)  =  V^i^<-^)     ^  exp  i  y   log  c^. -^  , 

valid  for  tl^  <  0  except  for  ^  =  0,  V)  >  -1.   These 
representations  make  it  evident  that: 


-8- 


y/i^)   Is  real  positive  on  the  positive  real  &x1b, 
C  =  4  =  0,  Y[    =  0, 

(^  +  1)  •V(^)  Is  real  positive  on  the  negative 
real  axis,  ^^  =  ^  <  0,  V]^  =  0. 

For  values  z  =  x  <  0,  y  =  0,  on  the  negative  real 

z-axis  we  may  take  as  the  path  (P   the  positive  real 

C-axls,  ^  =  ^  >  0,  ^n  =  0,  Prom  the  representations 

(7)  and  (8)  we  can  then  see  that  X.  (z)  and  hence  also 

2 

4  X^(z)  and  4  X^(z)  =  -  -^  X^(z)  are  real  for 

dz 
z  =  X  <  0,  y  =  0,  The  boundary  condition  II  is  therefore 

satisfied.   In  order  to  verify  that  the  boundary  condition 

^   is  satisfied  we  observe  that  the  integrand  in 

Is  regular  at  ^    =  -1.  Hence  we  can  deform  the  path  u 
into  the  negative  real  ^-axis,  C  =  ^  <  0,  Y|^  =  0,  if 
z  lies  on  the  positive  real  z-axis,  z  =  x  >  0,  y  =  0, 
Since  (  ^  +  1)  V(C)  is  positive  real  for  ^  =  ^  <  0, 
Y^  =  0,  we  see  that  (^  +  i)X  (z)  is  real  for  z  =  x  >  0, 
y  =  0,  The  same  statement  follows  similarly  for  7^  (z). 
The  boimdary  condition  I  is  therefore  satisfied. 

To  derive  the  behavior  of  X^  and  X^  at  infinity 

we  make  use  of  the  pole  of   V (C)  at  ^  =  -1.  Suppose  the 

i© 
point  z  =  re   lies  in  any  sector  excluding  the  positive 

real  axis,  -tt  <  0  <  e„„  <  0,   Then  ^  z  has  a  negative 

real  part  on  any  ray  C  =  pe   ,   X  =  constant  whose  angle 

\    lies  in  the  interval  -w<  ^  <   -5-  ©.«.".•  Since  any 

such  ray  in  the  lower  half -plane,  run  through  from  0  to  oo  , 

leaves  the  point  ^  =  -i  to  the  right.  It  is  admitted  as  a 

path  (P  •     If  we  now  let  z  approach  oo  it  is  clear  that 
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only  the  neighborhood  of  ^  =  0  will  contribute  asymp- 
totically to  the  integrals  (7)  and  (8).  By  (14)-  the 

1  ° 

-  w 

f\mctlon   V(^)  behaves  like  ^  at  the  origin; 

R  1 

consequently,   X  (z)  behaves  like  — •= —  for  large  values 

of  z  and  hence  vanishes  as  z  ->  oo  .   If ,  on  the  other 
1  Q 

hand,  z  =  re   lies  in  the  sector  -rr  <  ©„  <  ®  <  0»  we 

shift  our  path  to  any  ray  ^  =  pe   >  ^  =  constant,  with 

-^Tr-©„<A<-5#  and  then  argue  as  before.   In 

doing  so,  however,  we  cross  the  pole  at  C   -   -i«  The 

itr 

residue  contributed  by  the  pole  is  just  -r=.  e   e'   for 

iiT  J2 

T(^(z)  and  -L  ie"^  e"^^  for  X^(z)  as  seen  from  (12). 
^         /2 

Thus  the  formulas  (2)  and  (3)  are  justified. 

To  determine  the  behavior  of  the  functions   A(z) 

at  the  edge  z  =  0  we  must  establish  an  expansion  of 

V'(^)  for  large  values  of  ^  ,  Since 


.00 


dC. 


(16)  /  iog^„  -4^  =  ^-^(iogC+1)  -iC^dosC^)  +  ••• 
V  "  c  +1 


we  have  from  (14)^ 

00 


(17)   V(^)  =  C^   +|C"^(log^+l)  +  0(^-^  log^^) 


—2 
Thus   V  (^)  vanishes  at  infinity  like  ^        and  the 

integral 

(18)  v=.^fv(C)iC 

(p 
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is  absolutely  integrable  and  does  not  depend  on  which 
path  ^    Is  chosen.  Hence,  X^(z)  approaches  P  if  z 
approaches  0  through  values  for  which  the  path  (f    is 
admissible,  see  page  5.  Consequently,  if  the  fxmction 

X.^(2i)  is  defined  f or  z  =  0  by  X^(0)  =  P,  it  is  con- 
tinuous there. 

To  determine  the  value  of  P  we  consider  the  inte- 
gral 

(19)         '^  =  ^J  [cno  "  c(^^+i)"^]^^  » 

which  is  also  independent  of  the  choice  of  u   •  The 
term  subtracted  from   ^Ip^O    ^^   s°  chosen  that  the 

1    —2       y 

integrand  vanishes  at  infinity  like  —  ^    log  ^  ,  as 
seen  from  (17),  and  remains  integrable  at   C  =  0» 

If  one  takes  the  path  (F    in  (18)  and  (19)  as  the 
positive   ^  -8jcis,  it  becomes  clear  that  both  P  and  Q 
are  real. 

We  now  combine  the  two  integrals  to 

(20)   Q  +  iP  =^  I  (  (C+1)  V(C)  -  ^(^^+l)"^]d^  • 

Since  (C+i)  'V'(C)  ^^  regular  in  the  whole  lower  half -plane 
(^  ^  0)i  we  may  swing  the  path  of  integration  into  the 
negative    s  -axis,  taking  the  residue  due  to  the  second 
contribution  into  account.  We  obtain 


o     V 


(21)   Q  +  iP  = 


Now,  (^  +  i)  V(0  ^^  real  on  the  negative   ^-axis.  Hence 
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we  have 


(22)  ^  =  i  ' 


Q  =  ^  J     [(C+1)  no  -  c(^^+i)"^]  dc  . 


v-OO 

(23) 


Relation  (22)  yields  formula  (4)  since 

(24)  X^(z)  ->  I*   as   2  ->  0  , 

as  was  shown  above. 

We  now  apply  the  same  process  in  the  opposite 
direction.  We  consider  the  integral 

(25)  R  =^  J      I  C(r+  i)  HO 

o     >• 

-  [1  +  1  ^(^2  +  l)-^(log  t   +iTr  +  1)]  I 


dC 


which,  by  (17),  converges  absolutely  and  permits  swinging 
the  path  of  integration  back  to  the  positive  real  axis, 
except  for  a  residue  due  to  the  second  contribution.  We 
note  that  R  is  real  since  (C+  i)  V'(0  ^^^   1°S  4^  +  Itt 
are  real  on  the  negative  ^  -axis.  We  now  combine  R  with 
Q,  as  given  by  (25)  and  (23),  to 
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(26)  R .  iQ  =  ^  J     ((^:^  + 1)  v(C) 


-  [1  +  I  CiC^  +  i)"^(iog  C  +1)]  ]  dC 


and  swing  the  path  back  into  a  path  (T  •  Since 

(^^+1)  VCC)  ^s  regvilar  in  the  lower  half -plane,  we  obtain 

(27)  R  -  iQ  =  J^  I  {    if  +1)   V(C) 


^^n 


-[1+i  C(4^+l)"^(logr+l)]?d5  -|-^. 


The  integral  in  (27)  is  real  if  (T    is  taken  as  the 
positive   ^  -axis.  Consequently,  we  have 

(28)  Q  =  ^  . 

Prom  this  result  we  can  derive  formula  (5),  which 
gives  the  behavior  of  X  (z)  at  z  =  0,  We  first  observe 
that,  by  the  definition  (19)  of  Q,  relation  (28)  implies 


(29) 


^/^^'  [cw)-   C(C'-i)-"]d^  -^ 


if  z  approaches  0  through  values  for  which  u      is  admissible, 
Thus,  by  (8),  the  behavior  of  X  (z)  at  z  *=  0  is  reduced 
to  that  of 
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/  e^^  <C(C^  +  1)-1  dC     . 


In  order  to  discuss  this  integral  we  may  lead  the 
path  (F    first  along  the  positive   ^  -axis  from  0  to  any 
point  C   =   '^  f     '^  >  0,  and  then  to  infinity  in  such  a  way 
that  it  eventually  r\ins  along  a  straight  line  on  v/hich 
4z  is  real  and  negative.  Consequently,  we  may  write 

(30)  J  e^^iC^+ir^dl  =  J    e?^(C^+l)-^W 
j>  o 

'tOO 

+  1   e-V+z2)-Vd<r 

"-0(2 


-.<»< 


=  i       {e^^'DiC^+ir-CdC 


^ 


-  log  (-Z)  +^{l-e'"*^)    log  (z^-H^^z^) 


.00 

i  I    e"^log  (o-^+z^)d<r  . 


The  path  in  the  last  integral,  which  begins  at  the  point 
-oiz  in  the  upper  half -plane,  enters  infinity  along  a 
section  of  the  positive   <^-axis;   if  -o(z  lies  in  the  left 
half -plane  the  path  skirts  the  point  O"  =  -iz  so  as  to 
leave  it  to  the  right.   Letting  nov/  z  approach  0  we  obtain 

(31)   '  e^^(C^+l)"^  CdC  +  log  (-Z)  ->  /    e"^log  <Td<r    »  -C, 
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C  being  the  Euler  constant.  Consequently  we  have  from 
(8),  (29),  and  (31) 


(32)  A^(z)  -^  log  (-Z) 


1  «•  C 


as   z   ->  0,      Thus  formula   (5)    is   established, 

o 

Since  the  fimction  -X  (2)  is  the  derivative  of 
X  (z)*  the  formula  (5)  gives  at  the  same  time  a  refined 
behavior  of   X^(z)  at  z  =  0, 


(33) 


X^(z)  =  i  -  ^  {  2  log  (-Z)  -  (2-C)z  I  +  0(z), 


the  symbol  0(z)  representing  a  function  for  which 
z  0(z)  ->  0  as  z  -*•  0, 

Refined  asymptotic  representations  of  the  functions 
X  (z)  and   X  (2)  at  z  =  0  could  be  derived  by  continuing 
the  process  of  swinging  the  path  of  integration  of 
appropriate  integrals  forward  and  backward;  but  this  pro- 
cedure becomes  very  cumbersome.   Two  different  methods  of 
deriving  asymptotic  expansions  of  analytic  functions 
satisfying  the  two  boundary  conditions  I  and  II  on  the  two 
rays  of  the  real  ajcis  will  be  given  in  a  later  publication. 
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E:cplanation  of  the  figures 

The  figures  show  the  values  of  the  velocity  potential 
(pix,0)    =   Re  X(^)  along  the  dock,  x  <  0  and  along  the 
free  surface  x  >  0,  The  elevation  of  the  free  surface 

Is  given  by 

^  =  g"  (T  Cp(x,0)  slnfl-t,  X  >  0, 

the  pressure  under  the  dock  surface  by 

P  =  po-<p(x,0)  sino-t,  X  <  0. 

The   unit   of   length  is   chosen  as   g  <r"^. 

Figures    1   and  2  refer   to  the   v/ave  motion  described  by  the 

fTonction     X^(z)    and     X.^  {z)   respectively. 
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